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' We discuss two different types of issues concerning the quantization of Einstein- 

■ Rosen waves. First of all we study in detail the possibility of using the coherent 

j^j;^. states corresponding to the dynamics of the auxiliary, free Hamiltonian appearing 

J3 I in the description of the model to study the full dynamics of the system. For time 

periods of arbitrary length we show that this is only possible for states that are 
0^ ' close, in a precise mathematical sense, to the vacuum. We do this by comparing the 

quantum evolutions defined by the auxiliary and physical Hamiltonians on the class 
of coherent states. In the second part of the paper we study the structure of n-point 
Q^i functions. As we will show their detailed behavior differs from the one corresponding 

jl^ I to standard perturbative quantum field theories. We take this as a manifestation 

of the fact that the correct approximation scheme for physically interesting objects 
in these models does not lead to a power series expansion in the relevant coupling 
constant but to a more complicated asymptotic behavior. 

^ '■ PACS numbers: 04.60.Ds, 04.60.Kz, 04.62. +v 
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OO : I. INTRODUCTION 

o 

. ^ ^ Einstein- Rosen (ER) waves pLj] provide a very interesting toy model to discuss several 

^ ■ issues relevant for the quantization of general relativity jl, UH, d, [g], 0, Isl, M, [IdI, 11, [H, [3, 

. llJ, ll5|, [la, 113] ■ The reason behind this is the possibility of exactly describing the dynamics 



cr 



of the system both classically and quantum mechanically. This is true even after c oup ling 



some types of matter fields -massless scalars- to the gravitational degrees of freedom 
The main purpose of this paper is to discuss two different issues. The first is related to the 
problem of finding semiclassical states for the dynamics of the system, the second is to 
discuss the structure of important physical objects: the n-point functions of the model. 

A Hamiltonian description of ER-waves shows that the dynamics of this system is rather 
interesting owing to some unexpected features of the model. Probably the most striking 
one is that the Hamiltonian is bounded both above and below. This is a direct consequence 
of the fact that it can be written as a bounded function of the Hamiltonian for a free field 
theory [isl, l^- The origin of this free Hamiltonian can be traced back to the asymptotic 
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behavior chosen for the metric at infinity (in the 2+1 dimensional sense explained in jl8|). 
It is possible to gain a lot of information on the classical and quantum dynamics of the 
system by taking advantage of this functional dependence mentioned above. For instance, 
the quantization can be carried out by using a Fock Hilbert space defined by the auxiliary 
free Hamiltonian. This is so because the spectral theorem allows us to define the physical 
Hamiltonian once we construct the auxiliary one in a suitable Hilbert space. 

The first issue that we want to discuss here is related to the definition of semiclassical 
states for the full dynamics and the possibility of using the ones corresponding to the free 
auxiliary Hamiltonian as an approximate substitute. Coherent states play a very important 
role for systems of coupled harmonic oscillators -including free field theories that can be 
readily interpreted as models consisting of an infinite number of them. They display a very 
interesting behavior because they somehow bridge the gap between states with a purely 
quantum behavior and classical solutions to the equations of motion. For example, even 
though they display the characteristic dispersion of position and momentum observables, 
coherent states are minimal in the sense that the Heisenberg inequalities are saturated. 
Also the mean values of position and momenta evolve according to the classical equations 
of motion. For this reason coherent states can be considered as the best semiclassical states 
for linear systems. In the context of ER-waves, the large quantum gravity effects discovered 
by A. Ashtekar in ^ (see also [sl, 0, were analyzed by using coherent states. The main 
result in is that if one considers a coherent state for the quantum scalar field that 
describes the local degrees of freedom of an ER-wave and computes the relative uncertainties 
A$pO/(C)$(j of a certain relevant observable O -that can be interpreted as the quantum 
counterpart of a metric component of the ER-waves- one gets huge uncertainties even for 
coherent states with low (but no-zero) occupation number^ IIC'IP when C is peaked around 
a not too low value of the energy. Here, we want to discuss, from a dynamical point of view, 
the usefulness of coherent states as bona fide semiclassical states. As we will see, the class 
of coherent states can not be considered as semiclassical for the dynamics of the system 
and hence the behavior discussed in js], 0, 0, Isj is, in fact, rather natural from our point of 
view. The results of this paper should be considered as complementary to those presented 

in Baas. 

The problem of finding semiclassical states for systems different from the harmonic os- 
cillator is a very difficult one whose general solution is not known (in fact, even for such 
important systems as the hydrogen atom no such states have been found to date). As a 
consequence of this it is natural to expect that no coherent states -in the traditional sense- 
exist for the dynamics of the ER-waves (in fact this is a consequence of a simple exercise 
that can be carried out for the one-dimensional harmonic oscillator ^^). What we want 
to study here is to what extent the coherent states that do exist for the auxiliary dynamics 
can be used to derive meaningful information for quantum ER-waves. To this end, we will 
compare the states obtained by evolving a coherent state for the free Hamiltonian at a 
certain time to = both with the auxiliary dynamics t/o(t) and the full physical dynamics^ 
U(t). We do this by considering ||[/o(t)$c ~ U(t)^c\\, the projection of f/(t)$c on the 
coherent state $q labeled by the classical evolution of the classical initial data C. We do 
this both for small and large values of t. As we will see it is possible to quantify the periods 



^ Written in terms of the classical initial data C. 

^ Notice that we can use the same Hilbert space to describe the auxiliary and the full dynamics because 
the physical Hamiltonian is a function of a free Hamiltonian. 
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of time where free coherent states remain approximately coherent for the full dynamics. 
We will also show that for large values of t their distance becomes as large as possible (for 
orthogonal states). 

A second set of questions that we want to address concerns the relationship between 
the exact quantization of ER-waves and standard perturbative approaches. As the system 
can be exactly solved it is possible to devise efficient approximation schemes to extract 
physical information about it. We have done this in the past to discuss, for example, issues 



related to microcausality |10l . 1 121 . |13| . The main lesson that we have learnt from this type 
of analysis is that the asymptotic behavior of physically interesting objects is not captured 
by simple power series expansions in terms of the relevant coupling constant. Here we 
give an alternative way to understand this by looking at n-point functions. These are the 
building blocks used in standard perturbative analysis in QFT to obtain the S-matrix and 
discuss such important issues as renormalizability. As we will see the structure of n-point 
functions is such that one cannot expect a simple perturbative series to appear and, hence, 
the non-standard asymptotic behaviors found in previous works are natural in this setting. 

The paper is organized as follows. After this introduction we discuss in section [TTl the 
basic material needed to describe the canonical reduced phase space for the model. The next 
section IIIII is devoted to the construction of the one-particle Hilbert space for the auxiliary 
Hamiltonian. Section IIVI deals with the Fock quantization and the quantum dynamics of 
the system. We discuss several issues related to the definition of coherent states for the free 
auxiliary dynamics of the model in section |Vl We also study the asymptotics, both for small 
and large times, of some functions that measure the deviation of the states evolved with the 
full dynamics with respect to the states obtained by evolving the same quantum initial data 
with the free, auxiliary, dynamics. The main point discussed in section |Vl] is the structure 
of the ra-point functions. They can be exactly written in closed form. At every relevant 
order n we uncover a mixing property that characterizes their structure and shows that we 
are dealing with an interacting theory despite the fact that the main building block in its 
construction is a free model. We end with the conclusions and an appendix that summarizes 
some technical results concerning the asymptotic expansions that appear in the main body 
of the paper. 



II. CANONICAL REDUCED PHASE SPACE 

Einstein-Rosen waves describe vacuum solutions to the Einstein equations^ R^^^ = for 
a symmetry reduction of general relativity consisting of space times of the form (M^,gl^^) 
with four-dimensional metrics g^^^ having two hypersurface orthogonal, commuting, spatial 
Killing vector fields. The isometry group of these space-times is M x f/(l) and the metrics 
are regular at a symmetry axis. If we use a single global coordinate chart (x, y, z, t) on 
it is possible to take these Killing fields as {d/dzY and {d/daY '■= x{d/dyY — y{d/dxY- 
The cylindrical coordinates naturally associated to the previous cartesian coordinates allow 
us to write the four- dimensional metric in the form 

9ab = ^''"^[- e-^-(dt)„(dt), + (dr),(dr),] + r^e-'^ida)a{da)b + e'f'{dz)a{dzY 



^ Here and in the following we use Penrose's abstract index notation. 
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where — oo <t<oo, 0<r = \/ ^ y^, < cr < In, — oo < z < oo and the symmetry axis 
hes at r = 0. In the previous expression 7(t, r) = 7(t, a/x^ + y^) and r) = \/ + y"^) 
are smooth functions^ of {t,x,y) G M^, and 7oo(^) hm^^oo 7(^; Notice that e*^ is the 
norm of the translational Kilhng field (d/dz)^. The time coordinate t is chosen in such a 
way that the vector field {d/dtY is an asymptotic (for large r) unit Killing vector field for 
all the 2+1 dimensional metrics of the form 

9ah = e^{g'ab - e\dz)a{dz),) = -e^-^-(dt),(dt), + e\dr)a{dr), + r\da)a{da), (2.1) 

defined on the space of orbits of the translational Killing vector field. This choice makes sense 
0] because, under certain fall-off conditions for the fields 7 and 0, the Einstein equations 
imply that 700 is time independent and non-negative. The value of 

27r(l - e"^°°/^) G [0,27r) (2.2) 

represents the deficit angle of the asymptotically conical, 2+1 dimensional metrics, belonging 
to the class defined by (12. ip . We want to remark at this point that when 700 = this 
deficit angle vanishes and, hence, the asymptotic behavior is exactly given by the auxiliary 
Minkowskian background metric 

Vah = -(dt),(dt)b + (dr),(dr)b + r\<la)M<y)b ■ (2.3) 

This background metric will play an important role in the following. 

As we have mentioned above, the Einstein field equations R^^l = force 700 to be constant 
in t whereas the function 7 can be obtained in terms of (see, for example, jl, [l5|). In 
particular, once we fix some Cauchy data {Q, P) for at some initial time t = (say, 
to = 0), we can compute the quantity 

700 = ^^°° [P\r) + Q'\T))rdr , 

-that will be seen to be a constant of motion under the dynamics defined by the Hamiltonian 
given below- solve for 0(t, r) as a solution to the Einstein equations with initial data 

0(O,r)=g(r), e^°°/20(O,r) = P(r) r G [0, cx)) 

and, finally, obtain 

7{t,r) = - (e^-02(t,s) + 0'2(t,s))sds. 

As we can see, the local physical degrees of freedom of the ER waves are described by a 
scalar field 0. 

The dynamics induced by the Einstein equations on the field that we have just described 



^ These conditions guarantee the smoothness of the axiaUy symmetric scalar function at r = 0. In the 
following when we talk about smoothness in the axis we will refer to this condition. 
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admits a well-known Hamiltonian formulation [J] for which the Cauchy surfaces are the level 
surfaces of the asymptotic Minkowskian time coordinate t. In order to describe it let us first 
introduce C C C°^(]R^), the linear space whose points are smooth real functions on with 
rapid decay that depend on {x,y) G through r = ^ G [0, oo). The asymptotic 
conditions for the Cauchy data (5(r) and Pir) can be relaxed both at r = and r — ^ oo. This 
is important for the classical viewpoint but is irrelevant for the Fock quantization considered 
in this paper in the sense that the Fock space of quantum states turns out to be insensitive 
to the detailed choice of these asymptotic conditions. The canonical reduced phase space 
of the Einstein- Rosen waves T = ("P, tu), with points generically denoted as {Q,P) G T, is 
defined by endowing V = C x C with the standard (weakly) symplectic structure 

u{{Qi,P,),{Q2,P2)) (g2(r)Pi(r)-gi(r)P2(r))rdr. (2.4) 

The description of the classical dynamics in T is done in the form of an autonomous Hamil- 
tonian system {T,uj,h) with a Hamiltonian /i : T — > M that is defined in terms of the 
quadratic (free) auxiliary Hamiltonian^ 

ho{Q, P) := 7oo(Q, P) = \j^ (^'(0 + Q'\r)) rdr 

through a non-polynomial map 

h{Q, P) = 2 - 2 exp ( - i/io(g, P)) . (2.5) 

Notice that, in view of (12.21) . h{Q^ P) can be interpreted (up to a tt factor) as the deficit angle 
of the metrics (12.11) . The Hamilton equations derived from this non-quadratic Hamiltonian 
are 

Q = e-'^o«9,P)/2p 
P = e-^°(«'^)/2Ag 

where A denotes the Laplacian A : C ^ C acting on axially symmetric functions as 

F'ir') 

(AF)(r) ■=F"(r) + ^^ 

r 

Notice that the 1/r term in the previous expression originates in the axial symmetry of the 
system. 



^ We use units such that c = h = 8G3 = 1, where G3 denotes the effective Newton constant per unit length 
in the direction of the symmetry axis. 
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III. THE ONE-PARTICLE HILBERT SPACE ADAPTED TO THE 
ASYMPTOTIC STRUCTURE 

This section is devoted to the construction of the one-particle Hilbert space that we will 
later use to build the Fock space where the quantization of this system will take place. 
We start by pointing out that the Laplacian operator introduced above can be extended 
to a densely defined operator = —A on L^([0, oo), rdr). The operators G and a/O are 
self-adjoint and non-negative. As usual^ it can be employed to define a complex structure 
J : T — > T on the canonical phase space according to^ 



This complex structure is the restriction to the axisymmetric case of the standard complex 
structure adapted to the Poincare symmetry of the background metric (12.31) . It can be used 
to construct a complex vector space Tj whose points are exactly the same as the points of 
T and the multiplication by complex numbers^ x + iy E C is defined by 

{x + zy){Q,P) ■.= x{Q,P)+yJ{Q,P). 

It is possible to combine now u and J to define a positive definite sesquihnear form 

(■,-)j:TjxTj^C, 

((Qi, Pi), (Q2, P2))j = ^ouiJiQ,, Pi), (Q2, P2)) - ^^((Qi, Pi), (Q2, P2)) . 

providing us with a scalar product on Tj. The one-particle Hilbert space Tij of the ER 
waves is the Cauchy completion of (Tj, (■ , ■)j). 

There is another useful construction of the one-particle Hilbert space that uses Tc -the 
C- vector space obtained from T by considering complex functions in Cc C C°°([0, 00), C) 
with the standard multiplication by complex scalars- as the starting point. In order to see 
this we notice that the complex structure (13. ip can be diagonalized in T^. In fact, the 
vectors 

iC,TiVeC) eTc, CeCc, 

are eigenvectors of J corresponding to the eigenvalues ±i. Hence we can write Tc as the 
direct sum Tc = T_)_ © T_ where 

T± ■.= {{C,T^VeC)eTc\CeCc}. 
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Here we closely follow the ideas developed in [21 1 . 



^ In order to make sense of J it is necessary to restrict the domain of O > so that 1 / -y/O is well defined. 
To this end it suffices to consider functions F{x, y) such that their Fourier transform /(wi, ^2) vanishes in 
a neighborhood of zero. This guarantees that yju)l + w'^fiwi, ■W2) and f{wi, W2)/ y^wf + W2 are of rapid 
decay and smooth, even at (wi,W2) = (0,0), when f(wi,W2) is chosen smooth and of rapid decay. We 

will implicitly use this domain when needed. Notice however that there are many functions that can be 

2 2 

used to describe physical situations, for example the gaussian ~^ , that do not satisfy this restriction. 
^ Here x, y G M. As usual i — -\/^ is the imaginary unit. 
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It is clear that fl T_ = {0} and T_|_ = T_ where 

(c, := (c, ±ie^/^c). 

If we take a point (Q, P) G T there exists a unique C G Cc, given by 

such that 

{Q,P) = {C,-iVQC) + {C,-iVeC). (3.2) 

From equation (13. 2p it is easy to see that given the first component C = \{Q + iQ^^P) 
of (C, —i^/QC) G T+; the other can be then computed without any ambiguity. Hence, 
the one-particle Hilbert space Tij can be equally well described in terms of the complex 
functions C by using the following identification^ k : Tij ^ Ti 

C = k{Q,P):=]^(^Q + iQ-'^P^, (3.3) 
(Q, P) = K-^C = {C + C, -iVe{C -C)). (3.4) 

The map k is adapted to the complex form J in the sense that if k{Q, P) = C then 
Ko J(Q, P) = iC. According to the previous discussion we can build a Hilbert space 

n={C : ||Cf = < oo}, 

defined with the help of the scalar product 

(Ci,C2) := ]^uj{Jk-^C,,k-^C2)-'-uj{k-^C,,k-^C2), 

that is equivalent to the one-particle Hilbert space Tij. The complex structure in H. is 
diagonal {JC = iC for all C G 7i). At variance with the situation concerning the spaces T± 
(for which T+ = T_), if we work with H the operator is a conjugation in H, that is, Ti. 

is an antilinear map from Ti io Ti satisfying C = C. 

There are several mathematical structures that are easier to handle in Ti than in Tij. 
For example, the classical Hamiltonian that describes the dynamics of the Einstein-Rosen 
waves can be written now in terms of the scalar product and the operator G given above by 
noticing that 

hoin'^c) = \\Q-^cf = {c,Vec) 

h{K~^C) = 2-2exp(-(C, VeC)/2). 
It is interesting to point out here that the scalar product (■ ,■) onH can be written in terms 



^ At this point, H is just a linear space of complex functions C{r) that will become a Hilbert space once a 
scalar is introduced. 
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of the usual -^^^([0, oo),rdr) product 

POO 

Jo 

as (Ci,C2) = 2(Ci, v^C'2)l2- This allows to show that G is self-adjoint in Ti.. This can be 
seen by using the mode decomposition introduced at the end of subsection IIII Bl that allows 
us to write 9 as a multiplication operator. 



A. Classical dynamics in the one-particle Hilbert space 

Let us consider now the dynamics of the scalar field defined by the auxiliary free Hamil- 
tonian 

ho{Q, P) = \j^ (^'(^) + ^"(^)) = (*ith C = k{Q, P)) 

in the one-particle Hilbert space Ti. To this end we first write the free (linear) Hamilton 
equations in terms of the fields C, 

Ct = —iVOCt , with initial data Cq = C at t = 0. (3.5) 

The self-adjointness of 6 in allows us to write the general solution C° to the equation 
(13. 5p as the action of the unitary operator exp{—it^/&) on the initial data C, i.e. 

= exp{-itVQ)C . 

If = exp{—it^/&)Ca denotes the evolution from t = to an arbitrary time t of some 
initial data that we label as Ca, the unitary character of the classical evolution in Ti. implies 
that 

(CO,, = (Ci, C2) for alH e M . 

In particular the norms of the states remain constant in time. 

The non-linear dynamics defined by (12. 5p on the reduced phase space of the Einstein- 
Rosen waves can be also described in the one-particle Hilbert space Ti. The solution Ct at 
time t to the field equations with initial data C at t = can be written now as 

Ct = exp (^-zte-^o('^"''^)/2^) C . (3.6) 

Notice that Ct depends on the initial data C in a non-linear way. Hence, if we denote by 
Cit, C2t, and C(i+2)t the solutions corresponding to the initial data Ci, C2, and Ci + C2 it 
is clear that C(i+2)t 7^ Cu + C2t- From (13. 6p it is also evident that the evolution does not 
preserve the scalar product on 7i i.e. 

{Cu.C2t) ^ {Ci,C2) . (3.7) 
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However, it is easy to see from (13.61) that 

liail = foralHeM. 
This is not in conflict with (13.71) because Ct does not depend hnearly on C. 

B. Mode decomposition 

In the following we will find it convenient to work with a suitable mode decomposition 
adapted to the axial symmetry of our system. Given any F G C or Ce, we will use the 
following Fourier integral representation^'' 



1 

F{r) = J f{^)M^r) dw 



related to the two-dimensional Fourier transform 

{J^F){w) = {J^F){Vwl + wl ) = — ! F{^/x^ + y^) e-^(^'"i+2^"'2) (^^^^^^^ 

according to f{w) = y/2w{J^F){w). The functions f{w)/w, where w = ^/wf + belong 
to the class of C°°(]R^) in the two real variables {ui, 002) with rapid decay. Notice that with 
these conventions /(O) = 0. The action of the operator can be written in a nice way in 
this Fourier representation 

1 r°° 

{\/QF){r) = —= / wf{w)Jo{wr)dw. 
v2 Jo 

The mapping k connecting the {Q, P) and C descriptions of the one-particle Hilbert space 
also has a simple expression in this representation. In fact, given 

1 f°° 

Q{r) = J q{w)Jo{wr)dw 
1 f°° 

P{r) = —= / p{w)Jo{wr) dw 
v2 Jo 

then C = k{Q, P) (we will say that C and {Q, P) are K-related) if it can be written as 

1 f°° 

C{r) = —= / c{w)Jo{wr) dw , 



V2J0 

with 



, , 1 / , , ip(w) 



2 V w 
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Here Jo denotes the 0th order Bessel function of the first kind. 
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Finally the scalar product of the one-particle Hilbert space Ti. simplifies to 



r>00 

(Ci,C2) = / c^{w)c2iw) dw , \\Cf= I \c{w)\^dw. 
Jo Jo 

and also the free Hamiltonian can be written in the simple form 

POO 

hoin-^C) = / w\c{w)\'^dw. 



IV. FOCK QUANTIZATION AND QUANTUM DYNAMICS 

The one-particle Hilbert space H allows us to construct the Hilbert space J^si'H) used 
in standard approaches to the quantization of Einstein- Rosen waves as the symmetric Fock 
space 



J^^{n) = 07^®^'^ with := C, 



n=0 



where Ti^"'^ denotes the symmetrized tensor product of n copies of Ti. We will write the 
inner product in J-'s{'H) in the form (■ | ■). 

Following the rules of second quantization 2^ we extend certain operators from 7i to the 
symmetric Fock space jFg(7i). If we are given a unitary operator exp{iA) : Ti. —>■ Ti. written 
in terms of a self-adjoint operator A : T>{A) G Ti Ti, with domain 'D{A), we can promote 
them to the Fock space J-'s{'H). In particular, there exists a self-adjoint operator 

dT{A) : V{dT{A)) C J'siH) ^ J'siH) 

and a unitary operator r{iexp{A)) such that 

T{iexp{A)) = exp{idT{A)) : Ts{H) Ts{H) . 

The operator dr(y4) is called the second quantization of A and is defined by 



dr{A) := 0A(" 



n=0 

where 



:= 

and I denotes the identity operator on Ti. In particular we can use this procedure to 
construct the free auxiliary Hamiltonian. This is defined in the one-particle Hilbert space 
in terms of \/Q according to 

hoiQ, P) = ho{K-'C) = (C, V&C) . 
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Also, the classical evolution in "H is described in terms of the unitary operator exp(— zt-\/G)- 
Hence, the second quantization of exp(— ztA/B) and \/Q will give us the free quantum unitary 
evolution and the free quantum Hamiltonian. Explicitly, the auxiliary free Hamiltonian^^ 
Hq is defined on J-'s{'H) as 

Ho := dr(Ve) , 

and the quantum evolution operator is given by 

Uo{t) := r(^exp(-itv^)) = exp - itdT{VQ)^ = exp{-itHo) . (4.1) 

Notice that, given C®" G TY®"^" the action of the free Hamiltonian Hq can be read from the 
formula 

p " " 

(i/oC^«")(ri,...,r„) = — ^ / llc{w,)Jo{w,r,)dw,. 



i=l j=l 



Furthermore, notice that -ffolw = and hence, ii C E H C J-'s{Ti-) belongs to the domain 
of Hq, we have^^ 

POO 

{C\HoC)= w\c{w)\Mw = hoiC). 







In order to make sense of the quantum counterpart of the full (physical) classical Hamil- 
tonian 02.51) we make use of the quantum free Hamiltonian Hq and the functional relation 
between the free and physical classical Hamiltonians Hq and h. In particular 



h{Q, P) = E{hQ{Q, P)) := 2 - 2 cxp(-/io(g, P)/2) 
where E : [0, oo) — > [0, 2) is the function 

E{x) := 2 - 2 exp(-x/2) . (4.2) 



The spectral theorems then guarantee that the operator 



H := E{Hq) 

is a well defined self-adjoint operator on J^si'H). It is important to notice that H is not the 
second quantization of any self-adjoint operator on Ti. In particular, in spite of the fact that 
the restriction of H to the one-particle Hilbert space satisfies H\-)-i = E{\/Q), the quantum 
Hamiltonian H ^ dT(E(\/Q)). Hence the unitary operator evolution 



U{t) = exp{-itH) 

generated by H is not the second quantization of any unitary operator on the one-particle 
Hilbert space. This is not a surprise because, as we have discussed in section |Tll the full 



We use a lowercase ho to denote the Hamiltonian quadratic form in the one-particle Hilbert space and Hq 
for the free quantum Hamiltonian operator in the Fock space. 
In the following we use ho{C) to refer to hQ{K^^C). 
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classical dynamics is not even described by a linear operator in 7i. Finally, it is important 
to point out that the n-particle subspaces 7-^®"" of the Fock space are stable under the 
quantum evolution generated by H. At first sight this might seem striking because classical 
ER-waves are not stationary space-times and, in principle, one would expect particle creation 
effects. However, the asymptotic conditions j3] used to derive the Hamiltonian formulation 
discussed in section |Tl] restrict the class of ER-waves considered here to those metrics that 
are asymptotically Minkowskian in its 2 -|- 1 formulation (12. ip . In this context it is possible 
to use the preferred Fock quantization associated to the Minkowskian metric (I2.3P for which 
the particle creation effects are absent. 



V. COHERENT STATES 

As we have discussed in section [Til a vector in the one-particle Hilbert space C ETi can 
be thought of, through the identification fl3.3p - fl3.4p . as the Cauchy data {Q,P) at a given 
time for the scalar field that describes the degrees of freedom of an ER-wave. It is well 
known that there exists a family of quantum states £ ^si'H), parameterized by C G 7i, 
that behave semiclassically under the free auxiliary evolution. These are the coherent states 

oo 



n=0 ^ 



where C"^° = 1 G C and C"^" G 7^*^°" denotes the tensor product of n copies of the vector 
C E H. Notice that H^cll = 1 irrespectively of the value of ||C||. The scalar product 
of two coherent states and can be expressed in terms of the scalar product in the 
one-particle Hilbert space as 



($c, I = exp ( - Csf + 2lm(Ci, C2) 



in particular 



|(<l>cJ$c.)| = exp(^-^||Ci-C2f^ >0, for all Ci, C2 E H . 

The inner product {^Ci I ^02) never vanishes but |($Ci | ^C2)\ decreases when we increase 
the distance between the Cauchy data Ci and C2. The class of coherent states is closed 
under the free dynamics defined by fl4.ip . 

f/o(t)$c = exp{-ttHo)^c = ^cM-^tVe)c = ^c? ■ 

In other words, at any given time t the free quantum evolution of the coherent state as- 
sociated to the Cauchy data C is just the coherent state associated to the classical time 
evolution of these Cauchy data. 

For the full evolution the situation is, on the other hand, quite different because in this 
case, if C 7^ 0, the time evolution defined by the full physical Hamiltonian H is such that 

U{t)^c = exp{-itH)^c 7^ ^Ct = exp - ite-'^^^^^/^i^o) $c • 



13 



As we can see U{t)^c does not give the coherent state labeled by the classical solution Ct- 



Furthermore an argument similar to the one presented in [20[ for the harmonic oscillator 
shows that U(t)^c, with C 7^ 0, does not belong to the class of coherent states. The case 
C = is special because the coherent state $0 = 1©0©0©-- - G TsiTi) is both the 
Fock vacuum and the vacuum for the Hamiltonian H . It satisfies if$o = and hence 
U{t)% = $0. 

In the following we will give a quantitative measure of how the time evolution of the free 
coherent states deviates from the behavior that one would naturally demand for a bona fide 
coherent state. First, we will study the function 



Dcit) := \\Uomc-Um 



c\ 



|2 



2{<^c\ {1 - cos{t{E{Ho) - Ho)))<^c) 
00 ^ 

2 _ 2e-ll^ll' y -(C®" I cos (t(E(Ho) - iJo))C®") 



n=0 



that explicitly measures the distance between the states obtained by evolving a given coher- 
ent state with the free and the full dynamics. Second, we will consider the function 

Pc{t) := I Umc) = ($c I exp(zt(e-'^»(^)/2^o - i5(i^o)))$c) 

00 

l|2 



-\\c\\ 



?i=0 

that tells us how the full evolution of a coherent state defined by some Cauchy data deviates 
from the coherent state associated to the full classical evolution of the same initial data. 
This is done by studying the projection of one state onto the other. In particular we will 
consider the short and long time limits of Dc{t) and Pc{t)- The short time limit will give 
us information about how fast a coherent state of the auxiliary free dynamics ceases to be 
semiclassical. The large time limit will lend us some information about how far from each 
other these states are if we let them evolve for a sufficiently long time. 



A. Asymptotic behavior for short times 

The behavior of these functions for short times can be obtained from the following result 
that can be easily derived by using a Taylor expansion. 



Let C e P(e") C n, then 



k=0 

\k 



k=0 



{2k)\ 

k 



where, as usual, {0)q, = (\I'|C\I^) denotes the expectation value of the observable O in the 
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normalized state There are several cases that we have to analyze separately 

Dc{t) = e{{E{H,)-H,f)^^ + 0{t'). 
Re(PcW) = l-^((e-'^°(^)/2i/o-£^TOr)*c + 0(t'). 

First of all we see that -Dc;(0) = and -Pc'(O) = 1- Also, as expected, the short time 
asymptotic behavior is controlled by the energy E{Hq). For Dc{t) we see that, as long 
as we choose states $c such that {{E{Hq) — ifo)^)*c small the values of Dcit) will be 
approximately zero (they behave as a constant times t^). In an analogous way, those states 
giving a small value for {e~^°'^'-'^/'^Ho — E{Ho))^^ will force Pc{t) to remain close to one for 
a longer period of time. 

B. Asymptotic behavior for long times 

The study of the asymptotic behavior for i — > oo is not as straightforward as the previous 
one and requires some work. In this case we will use the stationary phase method to obtain 

the sought for asymptotic behaviors. In the following it will be useful to work with finite 
sums instead of infinite series so, for each A'" e N, we start by defining the truncations 

E{Ho)))C^^) . 

These functions Dc{t,N) and Pc(t,N) involve a finite number of terms and approximate 
the corresponding Dc{t) and Pc(t) uniformly^^ in t. This means that if we fix a certain 
element C gTC and s > there exists a natural number Nc{s) e N such that 

\Dc{t)-Dc{t,Nc{e))\<2e 
\Pc{t)-Pc{t,Nc{e))\<e 

irrespective of the value of t. This results allow us to work with the approximations given 
hy Dcit, N) and Pcit,N). 

Let us first consider the asymptotic behavior of the squared distance Dc{t). If C(r) is a 
continuous function given by the expression 

C{r) = —= \ c{w)jQ{wr)dw e H. 
v2 Jo 

the approximations provided by Dc{t, N) have the following asymptotic behavior for i — > oo 



Dc{t, N) := 2 - 2e-ll^ll' A(<^^" I cos {t{E{Ho) - 

n=0 

Pc(t,N) := e-ll^ll'y-^(C®"| exp(it(e-'^°(^)/2i/o- 

' ^ n.' 



n=0 



The auxiliary mathematical results presented in this section are proved in the appendix. 
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(see appendix 

/3e + l 

Dcit, AT) ~ 2 - 2e-ll^lP - B^e-^^^^T cos + 1)) (1) ^ , (5.1) 

where Be and /?c > 2 are real numbers depending on the chosen state C{r). As we can see 
Dc{t,N) approaches 2 — 2e~"'"" as l/^/fi (or faster). The distance remains small when 
— i> (i.e. when C is close, in the || ■ ||-norm, to the value G 7i that labels the Fock 
vacuum state $o) cind it approaches its maximum value^*^ when ||C|| — > oo (due to the 
exponential decay in ||C|p this limit is reached very fast). 

Finally let us consider Pc{t)- In this case, for a continuous C E H, the approximations 
Pc{t,N) have the following asymptotic behavior 

Pc{t, N) = e-"^"' + ^ exp{ttg{C))F{C, N) + 0{l/t) , 
vt 

where 

g{C) :=(/io(C)+2)e-'^°(^)/2_2 

and F{C, N) is a fixed factor that depends on C and A^. They can be obtained by applying 
the stationary phase method as explained in the appendix. The main conclusion that we 
draw from the asymptotic analysis that we have carried out in the large time limit t — >■ cx) is 
that the coherent states corresponding to the free dynamics do not behave as semiclassical 
states for the dynamics defined by the full Hamiltonian of the system as soon as ||C|| ~ 1. 
In particular the quantum evolution U(t)^c of the coherent state defined by the initial data 
C, with ||C|| ^ 1, and the coherent state $ct labeled by the classical evolution of C become 
almost orthogonal for large times. 



VI. QUANTUM FIELD OPERATORS AND N-POINT FUNCTIONS 

The main purpose of this section is to look at the problem of quantizing Einstein-Rosen 
waves and related models from a perturbative perspective. This is an interesting issue 
because we have an exact quantization in our hands and, hence, we can compare exact results 
with those obtained by suitable approximations. In fact, as we have discussed elsewhere 
0,1 Il2l . Il3[ |. if we use asymptotic methods to extract the physical behavior of the model in 



terms of the relevant coupling constant (related to the Planck length) we are led to behaviors 
that cannot be captured by the power series expansions that one expects to get from a 
perturbative approach. The ra-point functions play a very important role in quantum field 
theory. In fact, for the standard physical models, they are the key ingredients to construct 
relevant physical quantities such as the ^-matrix. We will try here to study the structure 
of the n-point functions and compare them with those obtained from familiar QFT's such 
as QED. As we will see the structure of these Green functions is not the standard one 
corresponding to interactions defined by field-dependent potential terms. This gives us a 
different perspective concerning the failure of standard perturbative treatments to deal with 
the types of QFT's considered in this paper. 



For unit orthogonal vectors the maximum valued of the norm of their difference is ^/2. 
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In the following we will use creation and annihilation operators a*{C) and a{C) respec- 
tively. These are labeled with vectors C G 7^ in the one-particle Hilbert space introduced in 
section IIIII They satisfy the usual commutation relations 

[a(C,),a*iC2)] = {C,,C2)ld, 

where Id denotes the identity operator on J^s(H). The conjugation : H ^ H introduced 
above allows us to define subspaces of 7i consisting of purely real or imaginary vectors 

nR:= {C en\C = C}, Hi := {C en\c = -C] . 

They are related by the equality 

Hi = iHr . 

Now, given / G Tin, we can define the field and momentum operators 0(/) and vr(/) in 
terms of annihilation and creation operators 

:= aq) + a*{f)=a{f) + a*{f) 
vr(/) := aitf)+a*itf) = -taif)+ta*{f). 

They satisfy the commutation relations 

[0(/i),vr(/2)]=2^(/i,/2)Id. 

It is possible to introduce a single operator T(C), labeled by C = ^{Q + iQ~2 P)^ to describe 
both the field and its canonically conjugate momentum 

T{C) := a(C) + a*{C) 

= \ {<Q) + a{iQ--2P) + a*{Q) + a^iiQ-'^P)^ 

= ^(0(g)+7r(0-ip)). 

The commutation relations for these operators are simply given in terms of the symplectic 
form ([MD by 

[T(Ci),T(C2)] = -2^(Ci,C2)Id. 

We study now the (Heisenberg image) time evolution of the T(C) from an initial instant 
of time t = to a generic time t both under the free auxiliary dynamics and the full dynamics 
introduced above. 

Free dynamics. In this case we can immediately see that 

T%t,C) := U,\t)TiC)Uo{t) = exp{itHo)TiC)expi-itHo) = T(C°) 

where C° = exp{—ity/Q)C is the free classical evolution of the Cauchy data defined by C. 
The fact that the free dynamics can be written in such simple terms reflects in the form of 
the n-point functions F° defined as the vacuum expectation values of products of T°(t, C) 
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for different instants of time and 

= (T(C°,JT(C2\)...T(C°,J)^,. 

In fact, it is well known that 

-^2n+l(^l) ^I'l ^2, C2', ■ ■ ■', hn+1, C'2„+l) = 

and F2n(ti,Ci;t2,C2', ■ ■ ■',t2n,C2n) can be written in terms of two-point functions 
F2{ti,Ci;tj,Cj). For example the four-point function is given by 

F^{ti,Ci;t2,C2,t2,C2;t4,C4) — F2°(ti, Ci; ts, (73)^2 (t2, C2; t4, C4) 

+ F°(ti,Ci;t4,C4)F°(t2,C2;t3,C3) 
+ -^2 (^1) ^ij ^2, C'2)-^2 (^3) C'a; ^4, C4) ■ 

Full dynamics. Let us discuss now the evolution defined by the full Hamiltonian of the 
system. In this case the (Heisenberg) time evolution of the operators T(C) is given by 

T{t,C) := U-\t)T{C)U{t) = exp(ii//)T(C) exp(-ii//) . 

The unitarity of the time evolution implies that the (equal time) commutation relations 
between the T(C) operators are independent of t, 

[T{t,C,),T{t,C2)] = -iu;{Ci,C2)ld. 

However, at variance with the free evolution, it is clear now that 

T{t,C) ^ T{Ct) = T(^exp ( - ite-''°^''~"^^/^Ve)c'^ 
because in this case the classical dynamics 

Ct = exp(-iie-'^°('^"^)/2v^)C 

has a non-linear depencence on the initial data C. This also means that it is not possible 
to find any Bogoliubov relation of the form 

U-\t)a{C)U{t) = a(AC) - a*(B^), for all t, 

for any pair of operators At and Bt defined on the one-particle Hilbert space. It is important 
to notice at this point that, despite the naive expectation, we have that 

U-\t)a(C)U{t) ^ a{exp{itE{^/e))C) . 
In fact the real situation is the following. The time evolution of the annihilation and creation 
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operators a*{C) and a{C) is given by expressions^^ of the form 

a{t,C) := U-\t)a{C)U{t) = exp{-ite-^°^^ E{Vef>'^')a(C) , (6.1) 
a*{t, C) := U-\t)a*{C)Uit) = a*{C) exp(zte-^"/' ® EiVef'^') . (6.2) 

This can be easily proved by using the identities 



[a{C),H'o'] = Hl^a{<c)2C) - H^a{C) 

or equivalently 

a(C)H^ = /7^a(etC') . 

In the previous formulas fl6.ip and fl6.2p we have used a notation that tries to convey the 
interplay between the C's that label the operators and the Hilbert space states upon which 
they act. Given a n-particle state V G 7i®°", with Fourier coefficients v{wi, . . . ,Wn), the 
vector 

exp(-ite^^»/2 ® E{VQf^^')a{C)V e 
has the following Fourier coefficients 

exp - ite-("'i+-+"'"-i)/^E(?i;))cHw(u7, wi, . . . , Wn-i) dw . 

This mixing, due to the interaction present in the system, introduces important complica- 
tions in the computation of n-point functions for n > 2 and makes it quite different from 
the free case. 

The n-point functions for the Einstein-Rosen waves considered here are defined as 

F„(ti, Ci; h, C2; . . . ; t„, Cn) := ($0 I T(ti, Ci)T{h, C2) ■ ■ ■ T(t„, C„)<l>o) . 

Owing to the fact that we are dealing with an effectively interacting model these ri-point 
functions behave very differently from the free case ones. The two-point function can be 
very easily computed in this case 

i^2(ti,Ci;t2,C2) = {%\ait,,C;)a*{t2,C2)%) = (a*(ti, Ci)$o I a*(t2, C2)$o) 
= (exp {ttiEiVe))C, I exp (^t2^(v^))C2) 



exp ( i{t2 — ti)E{w) ]ci{w)c2{w) dw 



These two-point functions have been studied in detail in [12|, |13[. Notice that when the 
coupling constant of the model^^ G = G^h is reintroduced {G = Gs in units ^ = 1) it 
appears in non-polynomial form trough the expressions involving the function E defined 
in (14.21) . explicitly E{w) = ^(1 — e'"^^^). This is a distinctive feature of the model that 
ultimately leads to behaviors that cannot be written as powers of the coupling constants. 



These have been derived in a sHghtly different form in [lo| . 

G3 is the effective Newton constant per unit length in the direction of the symmetry axis. 
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Likewise, the four-point function can be obtained in a direct way and is given by 

F4{tu Ci; C2; ta, C3; t4, C4) = (a*(t2, C2)a*(ti, Ci)$o I a*(t3, Cs)a*{t4, C4)$o) 

but, at variance with the situation for the free four-point function, F4 cannot be written as 
a sum of products of two-point functions. The difference hes in what we caU a mixing term 

{a*(t2, C2)a*(h, Ci)$o I a*(ts, Cs)a*(U, 04)^0) 

that cannot be written as a product of two-point functions. Exphcitly, 

{a*{t2, C2)a*{U, Ci)$o I a*{ts, Cs)a*{U, C4)$o) = 

^ L Uo ^^^^^^^'^ ^''P 0^^^ ~ t2e-'"'^^)E{w2) + ihe-'"^l''E{wi))) dw2 | x 

X Ci(wi) 03(101) exp ^ — iii£'(wi)^dwi -|- 



00 / /"OO 



02(^2)03(^2) exp - t2)e "'^/^E(w2)j dw2j x 

X Ci(wi)c4(u'i) exp ^^(^4 — ti)E{wi)^dwi . 

We will refer to this situation by saying that this last term has a mixing order of four. 
Notice again the non-trivial behavior of this function in terms of the coupling constant of 
the model. 

Finally it can be shown in general that F2n+i — whereas the 2n-point functions 

-^2n(^l) Ci] . . .] t2n, C2n) = 

= ($0 I a{ti,Ci){a{t2,C2) + a*(t2,C2)) ■ ■ ■ (a(t2n-l, ^2n-l) + a*(t2n-l,C2„-l))a*(t2n, C2n)$o) 

always have a term with maximal mixing order of 2n. In every case there are two extreme 
situations as far as the mixing order of the different terms is concerned. On one hand^^ we 
can have 

n n 

(^0 I W dihk-l, C'2fc-l)a*(^2fc, C'2fc)$o) = Y\. -^2(^2^-1, C2k-l', t2k, C2k) 
k=l k=l 

that can be written as a product of two-point functions. On the other hand 

n 2n n 2n 

{^olllaitk,^) H a*{ts,Cs)<^o) = {l[a*{tk,Ck)<^o\ H «*(^s,C,)$o) 

k=l s=n+l k=l s=n+l 

is maximally mixed (in fact the mixing order is higher than the maximum present for 2(n— 1)- 
point functions). 



^"^ With the aim of simphfying some expressions we will use the following notation for products of operators 
riLi^fc ■.= AiA2---An. 
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VII. CONCLUSIONS 



We have studied the quantization of Einstein-Rosen waves in the reduced phase space 
obtained by imposing the asymptotic flatness condition of 18|] and using an asymptotic, 
unit, timehke, Kilhng vector field to parameterize the time evolution. We have discussed 
two different types of issues that are relevant to understand quantum Einstein- Rosen waves. 
The first issue that we have considered is related to the semiclassical limit of the system. 
Specifically we have studied to what extent the coherent states corresponding to the free 
auxiliary dynamics of the model can be thought of as semiclassical states under the evolution 
defined by the full non-quadratic Hamiltonian. The conclusion that we draw from our 
analysis is that for short periods of time (with a length determined by the Hamiltonian 
Hq as expected on general grounds) the free coherent states can, indeed, be considered as 
semiclassical. In the long time limit we recover from a dynamical point of view the results 
of Ashtekar [5| about the existence of large quantum effects in the system. Specifically we 
see that the free coherent states with low occupation number (||C|| ~ 1) do not behave 
semiclassically and it gets worse and worse for larger values of ||C||. We have discussed this 
by studying the distance and the mutual projections of certain states obtained by considering 
the different time evolutions relevant in the model (the auxiliary one given by Hq and the 
full physical one). If we consider the squared distance, as a function of time, between the 
states obtained by evolving coherent states with the auxiliary evolution and the full evolution 
we see that it approaches the value 2 — 2 exp(— ||C||^). If ||C|| is very small this distance 
remains small whereas it becomes significant once ||C|| ~ 1. Finally for larger values of 
||C|| the distance approaches the maximum value for orthogonal states. Similar conclusions 
can be reached by studying suitable projections. We want to emphasize that we have not 
proved the impossibility of finding good semiclassical states for the model but only that not 
all the semiclassical states for the free auxiliary model can be considered as such for the full 
dynamics. An interesting open problem is to find a sufficiently large class of semiclassical 
states representing classical ER waves corresponding to arbitrary Cauchy data C. 

A second point that we have studied is the mathematical structure of the n-point func- 
tions. The main reason to do this is to get some information about the possible perturbative 
analysis of the system. We have seen that the structure of the n-point functions is different 
from the one corresponding to a free QFT. This is noteworthy because the formalism that 
we have used here relies on the fact that our model can be conveniently described in terms 
of an auxiliary free model. The structure of the ra-point functions, for which we are able 
to give closed form expressions, displays the distinctive features of an interacting model 
because they cannot be written simply in terms of two-point functions. Also the type of 
non-local interaction underlying the model shows up in the detailed form of these objects 
that differ from the ones obtained for familiar systems where the interaction is just given by 
a field-dependent potential. This is compatible with the known fact (discussed elsewhere) 
that the asymptotic approach to the study of physical observables for this model leads to 
expansions in terms of the relevant coupling constant (that can be interpreted as an effective 
Planck length) that are incompatible with any power series (isl . . 



We want to conclude by remarking that the present model can be exactly solved. In fact 
the exact evolution operator and their matrix elements can be exactly written. This means, 
in particular, that there is no need to separately consider the n-point functions to construct 
physical objects such as the S'-matrix. Of course n-point functions are interesting objects 
with important physical interpretations (see [lol, 15|) so it makes sense to understand how 
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they can be obtained as we have done here. 
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APPENDIX A: ASYMPTOTIC EXPANSIONS 

This appendix contains the proofs of several results used to obtain the asymptotic be- 
havior of Dcit) and Pcif) in the large t asymptotic limit studied in subsection IV B[ We will 
write them in the form of propositions. 

Proposition A.l. Given C E 7i and e > there is a number Nc{£) G N, independent of 
t G M, such that 



and 



\Dc{t)-Dcit,Nci6))\<2e 
\Pcit)-Pcit,Nc{e))\<e. 



The proof is based on the fact that it is always possible to find Nc{e) G N in such a way 
that 



\2n 



< e 



n=Ncie) 



and hence, 

\Dcit)-Dc{t,Nc{e))\ 



2e-"^ll' J2 I cos {t{E{Ho) - i/o))C®") 



n=Ncis) 



CO _ 

n=Nc{e) 
oo 

< 2e-"^il' y ■ II cos itiEiH^) - //o))C^"|| 

n=Nc{e) 

oo oo _ 

< 2e-ll^ll' y i^||c«-||2 = 2e""^ll'y i- 1| Cf "< 2^. 



n=Nc(£) 



n=Nc(e) 



22 



We also have 

\Pc{t)-Pc{t,Nc{e))\ = 



< e 



oo ^ 

n=Nc{e) 
oo ^ 

5Z 



n=Arc(£) 
oo 



(C®" I exp(it(e-^°(^)/2//o - £;(//o)))C®") 



n=Aro(£) " n=Nc(e) 

This result allows us to work with the truncations given by Dc{t, N) and Pc{t, N). 



Proposition A. 2. Let 



1 /"^ 

C(r) = / c{w)Jo{wr)dw G 7i , 
v2 Jo 



and let us assume \c{w)\'^ = Bcw'^'= + 0(w^'=+^) {w 0+), -Be 7^ and (3c > 2 real 
constants depending on the state C. Then the approximations Dc{t, N) satisfy, in the limit 

t — > OO; 



Dc{t, AT) ^ 2 - 2e-ll^ll^ - ^.e-H^H Y ( ^ j cos (^(/?. + 1)) (^^ 
>ls a consequence, Dc{t,N) approaches 2 — 2e~^^'^^^^ as or faster. 



4\ — 



c+l 



This result follows from 



^ 1 

Dc{t,N)^2-2e-\\^\\'J2-^^n, 

n=0 



where Aq :— 1 and 



= / COS ( t(y^ - -E'(y^Wj)) ) |c(wi)|^ • • • |c(wn)|^dwi • • -dWn 

i[0,oo)" V ^.^1 ^.^1 / 

-j^ POO rco 

= - / exp {it{z - E{z)))Gn{z, C)dz+- exp ( - it{z - E{z)))Gn{z, C) dz . 
In the last integral we have performed the following change of variables 

n 

{wi, Wn-1, Wn) ^ (wi, . . . , Wn-u z) , where 2: = , 
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and defined the functions Gn{z,C) as 

„ n—1 n—1 

Gn{z,C):= / \c{z -\2wj)\'^Y\\'^i^k)\'^^Wi- ■ -dWn-l. 

These befiave, in tlie z limit, as Gn{z,C) ~ ^j^j^ (3^ > 2 real 

constants that depend on the state that we have chosen to start with (notice that, as we 
pointed out in subsection [IT 



[Bl the Fourier coefficients c{w) satisfy c(0) = 0). Hence, the 
stationary phase method |23l] gives the following asymptotics for A^, n > 1, in the t ^ oo 
limit 

B^^ ^ n{(3, + 1) 



2 V 2 



cosQ(n(/3, + l))) {^^ 



We conclude that the main contribution to the asymptotic expansion of Dc{t,N) is given 
by 2 - 2e-ll^'l'(v4[5 + ) = 2 - 2e-ll^li'(l + A^). This way we finally get CT . 

Proposition A. 3. If C ETi (regular enough) then Pc{t,N) satisfies the identity 



Pc{t, N) = e-l'^"^ + ^.^^i^^p^FiC N) + Om {t oo) 

where 



g{C) := {ho{C) + 2)e-''°(^)/' - 2 
and F{C, N) is a function that depends only on C and G N. 



The proof is a straightforward application of the stationary phase method [23| to 



(C®" I exp(zt(e-'^«(^)/2/fo - ^(//o)))C^®") = / exp ( it{e-'''^'^^^'\ - E{z)) ) Gn{z, C) 



W°(^)/^G„(/io(C), C) '''^ ^'^^^^^ ^ ^^^^ + 0{t-^) , n> 1. 
Notice that 

(C®° I exp(2t(e-'^°(^)/2i/o - E{Ho)))C''^) = 1 . 
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